Dynamic Stabilization in The Double- Well Duffing Oscillator 
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Bifurcations associated with stability of the saddle fixed point of the Poincare map, arising from 
the unstable equilibrium point of the potential, are investigated in a forced Duffing oscillator with 
a double-well potential. One interesting behavior is the dynamic stabilization of the saddle fixed 
point. When the driving amplitude is increased through a threshold value, the saddle fixed point 
becomes stabilized via a pitchfork bifurcation. We note that this dynamic stabilization is similar 
to that of the inverted pendulum with a vertically oscillating suspension point. After the dynamic 
stabilization, the double-well Duffing oscillator behaves as the single-well Duffing oscillator, because 
the effect of the central potential barrier on the dynamics of the system becomes negligible. 
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A periodically driven double- well Duffing oscillator [Q , 
which has become a classic model for analysis of nonlin- 
ear phenomena, is investigated. It can be described in a 
normalized form by a second-order nonautonomous ordi- 
nary differential equations, 

• 3 = A cos cot, 



x + jx — x + x 



(1) 



where 7 is the damping coefficient, and A and u> are the 
amplitude and frequency of the external driving force, 
respectively. 

The Duffing equation (|lj) with negative linear stiffness 
describes the dynamics of a buckled beam |||| as well as 
a plasma oscillator |I| . Its regular and chaotic dynamics 
has been analyzed in great details by Holmes using both 
the theoretical techniques and the computer simulations 
The results of this work have been also confirmed 
in experiments by Moon || for a buckeled elastic beam. 
Since then, a number of authors have studied the forced 
double- well Duffing oscillator in the past two decades, 
and found rich dynamical behaviors Q such as the fractal 
basin boundary between coexisting competing attractors 
||, hopping cross- well chaotic state and so on. 

Here we are interested in the dynamical behaviors as- 
sociated with the saddle fixed point of the Poincare map, 
arising from the unstable equilibrium point of the poten- 
tial. One interesting behavior, associated with chaotic 
dynamics, is the homoclinic intersection of the stable 
and unstable manifolds of the saddle fixed point. In 
Ref. [^) , Holmes showed that as the forcing amplitude A 
is increased the stable and unstable manifolds intersect 
transversally, giving rise to homoclinic motions. How- 
ever, as A increases further, another interesting behav- 
ior, associated with stability of the saddle fixed point, 
occurs. When A passes through a threshold value, the 
saddle fixed point becomes stabilized via a pitchfork bi- 
furcation. We note that the best-known example of this 



dynamic stabilization is the inverted pendulum with a 
vertically oscillating suspension point H. To our knowl- 
edge, this is the first report on such a dynamic stabi- 
lization in the double-well Duffing oscillator. After the 
dynamic stabilization, the behaviors of the double-well 
Duffing oscillator closely resembles those of the single- 
well Duffing oscillator because the central barrier of 
the potential has no significant effect on the motion of 
the system. 

For the numerical calculations we transform the 
second-order ordinary differential equation (|l|) into a sys- 
tem of two first-order ordinary differential equations: 



y, y — — 7?/ + x — x + Acosujt. 



(2) 



These equations have a symmetry S, because the trans- 
formation 



-x, y 



-V,t^t+- [T(period) = — ], 

A UJ 



(3) 



leaves Eq. (||) invariant. If an orbit z(t)[= (x(y),y(t))] is 
invariant under S, it is called a symmetric orbit. Other- 
wise, it is called an asymmetric orbit and has its "conju- 
gate" orbit Sz(t). 

For the unforced case of A = 0, there exist a saddle 
equilibrium point at (x, y) = (0, 0) and a conjugate pair 
of stable equilibrium points at (x, y) = (±1, 0). However, 
as A is increased from 0, one symmetric saddle-type or- 
bit and two asymmetric attracting orbits with the same 
period 2tt/uj arise from the saddle equilibrium point and 
the two stable equilibrium points, respectively ||]. We 
note that they become the fixed points of the Poincare 
map P, generated by stroboscopically sampling the or- 
bit points with the external driving period T. Hereafter 
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we will denote the symmetirc saddle fixed point and the 
asymmetric stable fixed points by z* and z* , respectively. 
Here we are interested in the bifurcations associated with 
stability of the symmetric saddle fixed point. Its linear 
stability is determined from the eigenvalues, called the 
Floquct multipliers, of the linearized Poincare map DP, 
which can be obtained using the Floquet theory 

Since the Poincare map P with a symmetry S has a 
constant Jacobian determinant (det) less than unity (det 
DP = e~ lT ), the only possible bifurcations of periodic 
orbits are saddle- node (SN), pitchfork (PF), and period- 
doubling (PD) bifurcations pQ] . When a Floquet mul- 
tiplier passes through 1, a PF or SN bifurcation takes 
place. On the other hand, when it passes through —1, a 
PD bifurcation occurs. Hopf bifurcations are excluded. 

Each bifurcation curve in the parameter plane will be 
classified by a pair (p, q) invariant along the curve. Here q 
denotes the period of the orbit and p denotes the torsion 
number, which characterizes the average rotation number 
of the nearby orbits during the period q [jnj. The tor- 
sion number (normalized by the factor 2ir) at the PF or 
SN bifurcation curve becomes an integer. However, when 
crossing a PD bifurcation curve, not only the period but 
also the torsion number doubles from an odd multiple of 
1/2 to an odd integer. To keep p as an integer, we choose 
the pair of (p, q) for the PD bifurcation as that of the 
period-doubled orbit. We also note that the SN bifurca- 
tions may create symmetric and asymmetric orbits with 
the same (p, q) . To differentiate them, their symmetries 
are also used to classify the SN bifurcations. For the case 
of the symmetric (asymmetric) orbit, a letter "s" ("a") 
will be added in the third entry such as [p,q, "s" ("a")] 
to label the corresponding "symmetric" ("asymmetric") 
SN bifrucation curve. 

By varying the two parameters A and u>, we numeri- 
cally investigate the bifurcation behavior associated with 
stability of the symmetric saddle fixed point for a moder- 
ately damped case of 7 = 0.1. The associated bifurcation 
structure in the u> — A plane is shown in Fig. [I]. Note that 
the region, hatched with vertical lines, is just the stability 
region of the symmetric saddle fixed point. It is bounded 
by a lower PF bifurcation curve PF(0,1), denoted by a 
dashed curve, and by an upper symmetric SN bifurcation 
curve SN(l,l,s), denoted by a dotted curve. When cross- 
ing the lower PF bifurcation curve PF(0,1), the saddle 
fixed point becomes stabilized through a PF bifurcation 
by absorbing a pair of asymmetirc fixed points. As a re- 
sult of this dynamic stabilization, a symmetric stable or- 
bit with period 2tt/lo (fixed point for the Poincare map), 
encircling the unstable equilibrium point of the potential, 
appears. Note that this symmetric stabilized orbit cor- 
responds to the symmetric stable orbit, arising from the 
stable equilibrium point of the potential in the single- 
well Duffing oscillator. Hence the dynamical behavior 
after such a dynamical stabilization becomes essentially 
the same as that of the single- well Duffing oscillator ||. 
Such stabilized symmetric fixed point disappears at the 
upper symmetric SN bifurcation curve SN(l,l,s) by ab- 
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FIG. 1. Stability diagram of the symmetric saddle fixed 
point for (a) uj > oj r (— \/2) and (b) uj < tu r . The hatched 
region with vertical lines is just its stability region. The sym- 
bols SN, PF, and PD denote the saddle-node, pitchfork, and 
period-doubling bifurcation curves, respectively. Each curve 
is also labelled by a pair (p,q); p and q are the torsion num- 
ber and period, respectively. To differentiate the symmetric 
and asymmetric SN bifurcations, the letters "s" and "a" are 
also added in the third entry such as [p,q, "s"("a")]. For other 
details, see the text. 



sorbing a symmetric unstable fixed point, born at the 
lower symmetric SN bifurcation curve SN(l,l,s). 

For the unforced and undamped case, locally the Duff- 
ing oscillator near the two stable equilibrium points at 
(x, y) — (±1,0) behaves as a soft spring with the natu- 
ral frequency v2- Hence the main resonance occurs at 
u) — uJr(~ \/2) in the linear limit. For uj > co r the sym- 
metric saddle fixed point z* becomes stabilized at the 
PF bifurcation curve PF(0,1) in Fig. |l|(a) by absorbing a 
pair of asymmetric stable fixed points z*. However, for 
lu < u) T , the asymmetric stable fixed points z* disappear 
at the upper asymmetric SN bifurcation curve SN(l,l,a), 
denoted by a dash-dotted curve in Fig. |l|(b), through the 
collision with the asymmetric unstable fixed points born 
at the lower asymmetric SN bifurcation curve SN(l,l,a). 
After that, the symmetric saddle fixed point z* becomes 
stable at the PF bifurcation curve PF(0,1) by absorb- 
ing a pair of asymmetric stable fixed points born at the 
lower asymmetric SN bifurcation curve SN(l,l,a), which 
will be denoted by z* n . We also note that the upper and 
lower asymmetric SN bifurcation curves form a "horn" 
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FIG. 2. Bifurcation diagrams (plots of x vs. A) for (a) 
= 3.3 and (b) u — 3.0. Here the solid line denotes a stable 
fixed point, while the dashed line represent an unstable fixed 
point. The symbols denote the same as those in Fig. 1. The 
phase flow of the stabilized symmetric orbit for A — 10 is 
denoted by a solid curve in the inset in (a), and its Poincare 
map is represented by a solid circle. 

with a cusp at lo — Lo ri as in the asymmetric Toda oscil- 
lator ||. 

We now present the concrete examples of bifurcations 
associated with dynamic stabilization of the symmetric 
saddle fixed point. The bifurcation diagram and the 
phase-flow and Poincare-map plots are also given for clear 
presentations of the associated bifurcations. 

We first consider the case of lo > ui T . For to = 2ui T [~ 
2v2] , a subharmonic resonance occurs in which the asym- 
metric fixed points z* become unstable by a PD bifurca- 
tion. Note that the PD bifurcation curve PD(1,2), be- 
longing to the subharmonic resonance, becomes folded 
back at w = u>f(~ 3.23) [see Fig. [j] (a)]. Hence, when 
lo > u>i no PD bifurcations occur for the asymmetric 
fixed points z*. As an example, we consider the case 
of lo = 3.3. As shown in the bifurcation diagram in 
Fig. ^(a), the symmetric saddle fixed point z* becomes 
stable through a PF bifurcation for A — 9.255 • • • by ab- 
sorbing a pair of asymmetric fixed points z*. Then a 
stabilized symmetric orbit, encirling the unstable equi- 
librium point of the potential, appears, which is shown 
for A = 10 in the inset in Fig. ||(a). After this dynamic 
stabilization, the dynamical behavior becomes similar to 
that of the single-well Duffing oscillator, because the ef- 
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FIG. 3. Bifurcation diagram (plot of x vs. A) for ui = 1.2. 
Here the solid line denotes a stable fixed point, while the 
dashed line represent an unstable fixed point. The symbols 
also denote the same as those in Fig. 1. 

feet of the central potential barrier on the dynamics of the 
system becomes negligible. Such a stabilized symmetric 
fixed point also disappears for A = 18.105 • • • through a 
SN bifurcation by absorbing a symmetric unstable fixed 
point, born for A — 1.309 • • • via a symmetric SN bifur- 
cation. 

In the range of w r < to < Wf, the mechanism of dy- 
namic stabilization of the symmetric saddle fixed point 
is the same, except for the bifurcation behavior of the 
asymmetric fixed points z*. As an example, we consider 
the case of lo = 3.0. As shown in Fig. |2|(b), the asymmet- 
ric fixed points z* lose their stability for A = 0.748 • • • 
by forward PD bifurcations, but they become restabi- 
lized for A = 4.991 • • • by backward PD bifurcations. 
Note that the subsequent bifurcations, associated with 
dynamic stabilization, are the same as those for the above 
case of u> = 3.3. 

We next consider the case of lo < u T . An example for 
lo = 1.2 is shown in Fig. ||. Unlike the case of lo > ui T , the 
asymmetric fixed points z* disappear for A = 0.101 • • • 
by the asymmetric SN bifurcations, and then jump phe- 
nomena occur in which the small asymmetric fixed points 
z* are replaced by the relatively large asymmetric fixed 
points z* n , born for A ~ 0.053 75 via the asymmetric SN 
bifurcations. After that, the replaced asymmetric fixed 
points z* n play the same role for the dynamic stabilization 
as the small asymmetric fixed points z* do in the above 
case of lo > u) T . Hence the symmetric saddle fixed point 
z* becomes stable via a PF bifurcation for A = 1.550 • • • 
by absorbing z* n . 

Finally, we discuss the bifurcation behavior of the 
double-well Duffing oscillator after the dynamic stabi- 
lization of the symmetric saddle fixed point. We note 
that the stabilized symmetric orbit, encircling the unsta- 
ble equilibrium point of the potential, corresponds to 
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FIG. 4. Bifurcation diagram (plot of x vs. uj) for A — 10. 
Here the solid line denotes a stable fixed point, while the 
dashed line represent an unstable fixed point. The symbols 
also denote the same as those in Fig. 1. The phase portraits 
of the three large symmeric orbits, born by the symmetric SN 
bifurcations of type (l,l,s), (3,l,s) and (5,l,s), are shown for 
uj = 2.0, uj = 0.95598, and u = 0.57935, respectively. 

the symmetric stable orbit, arising from the stable equi- 
librium point in the single- well Duffing oscillator. Con- 
sequently, the double-well Duffing oscillator behaves as 
the single-well Duffing oscillator, because the central po- 
tential barrier has no significant effect on the motion of 
the system. As an example, we present a bifurcation 
diagram (plot of x vs. uj) in Fig. I| obtained by the 
frequency scanning for A = 10. The symmetric orbit 
stabilized through the PF bifurcation of type (0,1) dis- 
appears for uj = 2.646 • • • by a symmetric SN bifurcation 
of type (1,1, s), and then jumps to a large symmetric or- 
bit born for uj = 9.306 • • • via a symmetric SN bifurcation 
of type (l,l,s). A phase portrait of the large symmetric 
orbit for uj = 2.0 is shown in the inset. As uj is de- 
creased, large symmetric orbits with higher odd torsion 
numbers p (p = 3, 5, . . .), encircling the three equilibrium 
points of the potential, appear successively. The phase 
portraits of the large symmetric orbits with torsion num- 
bers p = 3, 5 are given for uj = 0.955 98 and 0.579 35 
in the insets, respectively. Note that the large symmet- 
ric orbits with higher torsion numbers have an increas- 
ing number of loops. Furthermore, each symmetric orbit 
with odd torsion number p loses its stability through a 
symmetry-breaking PF of type (p+1,1). Two such PF 
bifurcations of type (2, 1) and (4, 1) are also shown in 
Fig. 4. Then large asymmetric orbits with broken sym- 
metry undergo period-doubling cascades. We note that 
all these bifurcation behaviors are essentially the same as 
those in the single-well Duffing oscillator (refer to Fig. 1 



in Ref. 0). 

To confirm the above numerical results, we also con- 
structed the electronic analog simulator of Eq. (|l|) with 
the usual operational amplifiers and multipliers, and 
made an analog study. Dynamic stabilization of the sad- 
dle fixed point was thus observed experimentally. The 
details on the experimental results will be given elsewhere 
PI- 
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